FERMILAB-PUB- 1 1- 163-PPD, BNL-94488-2010- J A, arXiv:hep-ex/l 104.0344 



First Direct Observation of Muon Antineutrino Disappearance 

P. AdamsonJ C. Andreopoulos.^o D. J. Auty,^^ D. S. Ayres,i C. Backhouse/^ G. Barr,i* M. Bishai,^ A. Blake,^ 
G. J. Bock/ D. J. Bochnlcin/ D. Bogert/ S. Gavanaugh,^ D. Gherdack,^^ S. Ghildress/ B. G. Ghoudhary/ 
J. A. B. Goelho,^ S. J. Goleman,^^ L. Gorwin/^ D. Gronin-Hennessy,!^ I. Z. Danko,!^ J. K. de Jong,!^ 
^ N. E. Devenish,24 M. V. Diwan,^ M. Dorman,!^ G. O. Escobar,^ J. J. Evans," E. Falk,24 Q. J. Feldniaii,^ 

(3 ! M. V. Frohne/° H. R. Gallaglier,^^ R. A. Gonies,^ M. G. Goodman/ P. Gouffon/i N. Graf," R. Gran,!^ N. Grant,^^ 
(N ; K. Grzelak,28 A. Habig,!^ D. Harris,^ J. Hartnell,^^. 20 Hatcher,^ A. Himmel,'* A. Holin,!^ G. Howcroft,^ 

tJo' X. Huang,! J. Hylen,^ J. Ilic,^^ G. M. Irwin,23 Z. Isvan,!^ D. E. Jaffe,^ G. James,^ D. Jensen,^ T. Kafka,^^ 

^3 ; S. M. S. Kasahara,i5 G. Koizumi,^ S. Kopp,^^ M. Kordosky,^^ A. Kreymer,^ K. Lang,^^ G. Lefeuvre,^* J. Ling,3.22 
<^ ■ P. J. Litchfield,i5'20 L. Loiacono,26 P. Lucas, ^ W. A. Mann,^^ M. L. Marshak,!^ Mayer.^^ A. M. McGowan,i 
'. R. Mehdiyev,26 J. R. Meier,!^ M. D. Messier,!^ W. H. Miller,!^ S. R. Mishra,^^ J. Mitchell,^ G. D. Moore,^ 

J. Morfm,^ L. Mualem,^ S. Mufson,!^ J. Musser,!^ D. Naples,!^ J. K. Nelson,^^ H. B. Newman,"* R. J. Nichol,!^ 
T. G. Nicholls,20 J. A. Nowak,!^ J. p. Qchoa-Ricoux,* W. P. Oliver,^^ M. Orchanian,'* R. Ospanov,^^ J. Paley,^'^^ 
X'; R. B. Patterson,-* G. Pawloski,23 G. F. Pearce,^^ D. A. Petyt,!^ S. Phan-Budd,i R. K. Plunkett,^ X. Qiu,^^ 

^ ; J. Ratchford,26 T. M. Raufer,^^ B. Rebel,^ P. A. Rodrigues,!^ G. Rosenfeld,22 H. A. Rubin," M. G. Sanchez,". L 9 

J. Schneps,^'^ P. Schreiner,i P. Shanahan,'^ A. Sousa,^ P. Stamoulis,^ M. Strait,*'^ N. Tagg,*'^ R. L. Talaga,i 
^ ■ E. Tetteh-Lartey,25 J. Thomas,** M. A. Thomson,^ G. Tinti,*® R. Toner,^ G. Tzanakos,^ J. Urhenn,*^ p. Vahle,^^ 
B. Vhen,3 A. Weber,*^ R. G. Webb,^^ G. White,** L. Whitehead,^ S. G. Wojcicki,23 T. Yang,23 and R. Zwaska^ 

(MINOS GoUaboration) 

^ , ^ Argonne National Laboratory, Argonne, Illinois 60439, USA 

' ^Department of Physics, University of Athens, GR-15771 Athens, Greece 

' ^ Brookhaven National Laboratory, Upton, New York 11973, USA 

CO , ^ Lauritsen Laboratory, California Institute of Technology, Pasadena, California 91125, USA 

' ^Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 OHE, United Kingdom 

^ \ ^ Universidade Estadual de Campinas, IFGW-UNICAMP, CP 6165, 13083-970, Campinas, SP, Brazil 

, ''Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 

^Instituto de Fisica, Universidade Federal de Goids, CP 131, 74001-970, Goidnia, GO, Brazil 
^Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 
"'Holy Cross College, Notre Dame, Indiana 46556, USA 

_^ I ^'Physics Division, Illinois Institute of Technology, Chicago, Illinois 60616, USA 

' '^Indiana University, Bloomington, Indiana 47405, USA 

}_( , ' 'Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011 USA 

■ ''''Department of Physics and Astronomy, University College London, Cower Street, London WCIE 6BT, United Kingdom 

'^University of Minnesota, Minneapolis, Minnesota 55455, USA 
'^Department of Physics, University of Minnesota - Duluth, Duluth, Minnesota 55812, USA 
'''Otterbein College, Westerville, Ohio 43081, USA 
'^ Subdepartment of Particle Physics, University of Oxford, Oxford 0X1 3RH, United Kingdom 
'^Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
^"Rutherford Appleton Laboratory, Science and Technologies Facilities Council, 0X11 OQX, United Kingdom 
^' Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 05315-970, Sao Paulo, SP, Brazil 
'"^'"^ Department of Physics and Astronomy, University of South Carolina, Columbia, South Carolina 29208, USA 

'^'Department of Physics, Stanford University, Stanford, California 94305, USA 
'"^^ Department of Physics and Astronomy, University of Sussex, Falmer, Brighton BNl 9QH, United Kingdom 

Physics Department, Texas A&M University, College Station, Texas 77843, USA 
Department of Physics, University of Texas at Austin, 1 University Station C1600, Austin, Texas 78712, USA 
^''Physics Department, Tufts University, Medford, Massachusetts 02155, USA 
Department of Physics, Warsaw University, Hoza 69, PL-00-681 Warsaw, Poland 
Department of Physics, College of William & Mary, Williamsburg, Virginia 23187, USA 

This Letter reports the first direct observation of muon antineutrino disappearance. The ML 
NOS experiment has taken data with an accelerator beam optimized for 77^ production, accu- 
mulating an exposure of 1.71 x 10^*^ protons on target. In the Far Detector, 97 charged cur- 
rent 77^ events are observed. The no-oscillation hypothesis predicts 156 events and is excluded 
at 6. So-. The best fit to osciUation yields |Am^| = [3.36l{];4[5(stat) ± 0.06(syst)] x 10"^ eV^ 
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sin^(2^) — 0.861q (stat) ± O.Ol(syst). The MINOS and measurements are consistent at 
the 2.0% confidence level, assuming identical underlying oscillation parameters. 

PACS numbers: 14.60.Lm, 14.60.Pq, 14.60.St 



Observations by many experiments provide compelling 
evidence for neutrino oscillation This oscillation, 

a consequence of the quantum mechanical mixing of the 
neutrino mass and weak flavor eigenstates, is governed by 
the elements of the Pontecorvo-Maki-Nakagawa-Sakata 
matrix [lo| . parameterized by three mixing angles and 
a CP phase, and by two independent neutrino mass- 
squared differences. As the measurement precision on 
oscillation parameters improves, so does the potential for 
observing new phenomena. In particular, measured dif- 
ferences between the neutrino and antineutrino oscilla- 
tion parameters would indicate new physics. CPT sym- 
metry, one of the most fundamental assumptions under- 
lying the standard model, constrains the allowed differ- 
ences in the properties of a particle and its antiparticle 
111 ] and requires that their masses be identical. This 
symmetry has been extensively tested in other sectors, 
most notably the kaon sector [l^l. Additionally, neu- 
trinos passing through matter could experience nonstan- 
dard interactions [13| that alter the and disappear- 
ance probabilities and, thus, the inferred oscillation pa- 
rameters [l^ . 

The MINOS experiment has used a beam to mea- 
sure the larger (atmospheric) mass-squared difference 
|Am2| (2.32^^:03) X 10-3 eV^ and the mixing angle 
sm'^{2e) > 0.90 (90% confidence limit [C.L.]) through 
observation of disappearance 0, [isj . The correspond- 
ing antineutrino oscillation parameters are much less pre- 
cisely known. 

This Letter describes the first direct observation of 
disappearance consistent with oscillation, yielding the 
most precise measurement to date of the larger antineu- 
trino mass-squared difference. The only previous mea- 
surements from 17^-tagged samples, providing weak con- 
straints, come from the MINOS atmospheric neutrino 
sample [lij and an analysis of the 77^ component of the 
MINOS data sample [H IH. The strongest indi- 
rect constraints come from a global fit dominated 
by Super-Kamiokande data which measure the sum of 
atmospheric and interaction rates. 

For this measurement the NuMI beam line [2^ was 
configured to produce a 77^j-enhanced beam. The current 
in the magnetic horns was configured to focus negative 
pious and kaons produced by 120 GeV protons incident 
on a graphite target. Most mesons travel along a 675 m 
long decay pipe, filled with helium at 0.9 atm, and decay 
to produce a F^-enhanced beam with a peak energy of 
3 GeV (see Fig. [J). Interactions of z^^ comprise a frac- 
tion of all charged current (CC) events in the MINOS 
detectors which rises from about 21% below 6 GeV up 
to about 81% at 20 GeV, in the case of no oscillation. 



The data set in this paper corresponds to an exposure of 
1.71 X 10^° protons on target (POT). 

The MINOS experiment uses two similar detectors lo- 
cated 1.04km [Near Detector (ND)] and 735km [Far De- 
tector (FD)] from the NuMI target. The CC interac- 
tion rate as a function of reconstructed energy is mea- 
sured in each detector. The measured FD energy spec- 
trum is compared to that predicted using the ND data. 
In this comparison, many sources of systematic uncer- 
tainty cancel. Antineutrino oscillation causes a deficit in 
the FD with an energy dependence, in the approximation 
of two-flavor mixing, of 



= 1 -sin^(26 



^ i 1.267 ArrfiL\ 
I sm- I 1 (1) 
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where L [km] is the distance from the point of antineu- 
trino production, E [GeV] the 17^ energy, Am^ [^V^] the 
antineutrino mass-squared difference and 9 the antineu- 
trino mixing angle. 

The MINOS detectors [2l[ are tracking calorimeters, 
formed of planes of steel interleaved with planes of scintil- 
lator. The scintillator is divided into strips with a width 
of 4. 1cm. In CC interactions, P'p(t'p)+A^ /^+(/i^)+A', 
a hadronic shower (X) and a muon track may be ob- 
served. The hadronic energy is measured by summing 
the amount of light produced in the scintillator. Muon 
energy is measured by range for contained tracks or, for 
exiting tracks, by the curvature in a 1.4 T toroidal 
magnetic field. The incoming neutrino energy is recon- 
structed as the sum of the hadronic and muon energies. 
For the data presented in this Letter, the fields in both 
detectors focus /i+ and defocus fi~ , allowing the separa- 
tion of I7p and Ufj, CC interactions on an event-by-event 
basis. 

A sample of 77^ CC interactions is isolated by identi- 
fying the presence of a positively charged track. Neu- 
tral current (NC) interactions produce only a hadronic 
shower at the vertex. Similarly, CC interactions of Ve 
and Ve (which correspond to 2.0% of all CC interac- 
tions at the ND) produce only showerlike activity. The 
main background arises from tracks reconstructed out 
of shower activity. This background is reduced 0, [2^ 
by a method which uses four variables to identify the 
presence of an isolated track with muonlike energy de- 
position. These four variables are the track length, the 
average pulse height per plane along the track, the trans- 
verse energy deposition profile of the track, and the fluc- 
tuation of the energy deposited in scintillator strips along 
the track, and are combined using a fc-nearest-neighbor 
(fcNN) algorithm [2^ to produce a single output vari- 
able. The position of the selection cut on this variable 
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FIG. 1: The reconstructed energy spectra of events in the 
Near Detector classified as charged current interactions, sep- 
arated according to the reconstructed charge of the track. The 
events with a negatively charged track are not used in the os- 
cillation analysis. The shaded bands represent the systematic 
uncertainty on the simulation. 



is tuned to optimize the statistical sensitivity to |Am^|, 
yielding the same selection criterion as for the MINOS 
analysis [H. 

The charge of reconstructed muon tracks is determined 
by analyzing the curvature of the track in the magnetic 
field \M. Fi gurc [1] shows the reconstructed energy of se- 
lected CC events in the ND, separated according to the 
measured track charge sign. The events reconstructed 
with a negatively charged track consist primarily of z^^ 
CC interactions, and are removed from further analysis. 
Events with a positively charged track form the selected 
CC sample, and are used to predict the expected en- 
ergy spectrum at the FD. Below 6 GeV, where the ma- 
jority of the oscillation signal is expected, the selected 
CC sample at the ND has a purity, obtained from 
the simulation, of 98% (the background consisting of 1% 
NC events, 1% i/fj^ CC events). Above 6 GeV the pu- 
rity is 88%, and the contamination is primarily CC 
events; higher momentum muons follow a less curved 
path, giving a greater probability of charge misidenti- 
fication. The total CC reconstruction and selection 
efficiency is 93%. 

The measured ND energy spectrum is used to predict 
the FD spectrum, as previously [H d, [l^. This proce- 
dure is particularly effective in mitigating sources of sys- 
tematic uncertainty which affect both detectors similarly. 
For example, uncertainties on the neutrino flux and cross 
sections dominate the systematic error band on the ND 



FIG. 2: The distribution of the sign of the reconstructed 
charge divided by the momentum of selected muon tracks 
in the Far Detector. The simulated distribution is shown in 
the case of no oscillation and oscillation assuming the best 
fit parameters from Ref. 0] and parameters from this 
analysis. 



energy spectrum, shown in Fig. [1] but have a negligible 
impact on the oscillation measurement. 

The production of hadrons in the NuMI target is con- 
strained by fits to the ND data [l|, Q . These fits use data 
from the i/^ beam to determine the tt and K yields as 
a function of their transverse and longitudinal momenta 
at production. Recent measurements [25| of the ratio of 
7r+/7r~ yields are included as constraints in these fits. 
This tuning procedure improves agreement between the 
simulated ND energy spectrum and the data, but does 
not significantly affect the predicted FD energy spec- 
trum. Uncertainties on the modeling of the beam have 
a negligible effect on the predicted FD energy spectrum, 
and are accounted for in the oscillation measurement. 

The same event selection criteria are used in both de- 
tectors. The FD data selection was determined using 
simulation and ND data, before the FD data was exam- 
ined. All FD events passing the fcNN selection are shown 
in Fig. [21 distributed according to the sign of the recon- 
structed track charge, divided by the momentum. The 
figure shows good modeling of track charge identification. 
Events with a negatively charged track are minimally af- 
fected by oscillation due to their higher mean energy, and 
are removed from further analysis. 

The systematic uncertainty on the oscillation param- 
eters is much smaller than the statistical uncertainty. 
The sources of systematic uncertainty are very similar 
to those discussed for the MINOS analysis Q. An 
additional uncertainty is estimated on the level of 
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FIG. 3: Comparison of the measured Far Detector 77^ CC 
energy spectrum to the expectation in three cases: in the 
absence of oscillation, using the oscillation parameters which 
best fit this 77^ data (for this case, the total expected back- 
ground is also indicated) , and using the best-fit oscillation 
parameters measured by MINOS [3]. 



CC background in the selected 17^ CC event sample; be- 
low 6 GeV, the purity of the selected 77^ CC sample is 
known to within 1%. To evaluate this uncertainty, the 
simulated background is scaled until the total number 
of simulated and data events match in the background- 
enhanced set of events which fail the fcNN selection cri- 
terion. This scale factor is taken as the uncertainty on 
the level of background in the selected 17^ CC sample. 
The total systematic uncertainty on the measurement of 
|Am^| is +0.063 — 0.060 x 10~'^eV^; on the measurement 
of sin^(26') the total systematic uncertainty is ±0.012. 

Using the prediction obtained from the ND data, 156 
selected 77^ CC events with energy below 50 GeV are ex- 
pected in the FD in the absence of oscillation while 97 
events are observed. The energy spectra of these FD 
events are shown in Fig. [3l A clear energy dependent 
deficit is observed, showing the first conclusive evidence 
for 17^ disappearance consistent with oscillation in a 17^- 
tagged sample. The no-oscillation hypothesis is disfa- 
vored at 6.3 standard deviations. 

Oscillation is incorporated into the predicted energy 
spectrum according to Eq. [TJ Comparing the prediction 
to the data using a binned log likelihood, oscillation pa- 
rameters are found which maximize the likelihood. These 
f\ = [3.36t[!;4j5(stat)±0.06(syst)] xlO-^eV^ and 



are 



sin^(26l) = 0.86To:i2(stat)±0.01(syst)jand are consistent 
with all other previous direct limits |16l4l8l |. The predic- 
tion for oscillation with these best fit values is shown in 
Fig. O and includes 2 NC events, 5 CC events and 0.3 



FIG. 4: Allowed regions for the 77,j oscillation parameters from 
a fit to the data in Fig. [3] including all sources of sy stematic 
uncertainty. Indirect limits prior to this work [ly] and the 
MINOS allowed region for ly^^ oscillation 0| are also shown. 



Ur CC events. 

The confidence limits on the oscillation parameters, 
shown in Fig. 21 are calculated using the Feldman- 
Cousins technique [l^ , in which the effect of all sources 



of systematic uncertainty is included [27|, |28[ . Figure |4] 
compares these limits to the previous best limit [l9| . 

MINOS has measured the z^^ oscillation parameters 
to obtain a best fit of [Am^l 2.32 x lO^^eV^, 
sin^(26') = 1.0 0. Assuming that muon antineutrinos os- 
cillate with these parameters, 110 selected events are ex- 
pected in the FD below 50 GeV. This expected energy 
spectrum is shown in Fig. |31 denoted as 'V^ best fit" . 

Figure [3] compares the MINOS measurements of the 
ly^ and 17^ oscillation parameters. In both measurements 
the purity of the event samples in the oscillation region 
is high. Below 6 GeV, there is no more than 3% t'^ 
CC contamination in the 17^ CC sample and vice-versa. 
Therefore the measurements of the and j7^ oscillation 
parameters are nearly independent. Since the 17^ mea- 
surement is heavily statistically limited, the impact of 
correlated systematic uncertainties is negligible. 

In a joint fit to the data used in the MINOS i/^ and 
Vfj, measurements, assuming identical and 77^ oscil- 
lation parameters, the best fit parameters are |Ar77^| = 
lAm^l = 2.41 X 10-3 eV^ sin2(2^) = sin2(26') = 0.97. 
The significance of the difference in likelihood between 
this best fit and the fits to the individual i/^ and 17^ data 



sets is evaluated using a Feldman-Cousins approach [27 1. 
Ten thousand simulated experiments are generated as- 
suming the joint best fit oscillation parameters above, 
and include all sources of systematic uncertainty. The 
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difference in likeliliood between tlie joint and individ- 
ual and fits is recorded for eacli experiment, and 
the fraction of simulated experiments with a difference 
in likelihood larger than that observed in the data is a 
measure of the significance of the observed difference. 
Assuming identical and oscillation parameters, the 
probability that the MINOS measurements of the and 
Vp, parameters would be more discrepant than those ob- 
served is 2.0%. 

A thorough search for sources of mismodeling that 
could provide a difference between the and mea- 
surements was performed. The only change between 
and Vp running modes was the reversal of the directions 
of the current in the focusing horns of the beam and of the 
magnetic fields in the detectors. Very similar data anal- 
ysis procedures are used in both modes, with the same 
reconstruction code and similar selection criteria used to 
obtain the event samples, and the same technique used 
to obtain the FD predictions. These similarities make 
the comparison of and measurements robust and 
limit the possible sources which could generate a spurious 
difference. 

The and analyses differ in that the F^-enhanced 
beam contains a significant component (which is ef- 
fectively removed by the selection cuts) . Figure [1] shows 
that this component is at high energy, away from the os- 
cillation signal region, and therefore any residual contam- 
ination has little effect on the oscillation measurement. 
Figures [T] and [2] show the z/^ CC component to be well 
modeled in both detectors in both shape and normaliza- 
tion. All FD events were scanned by eye to ensure the 
selection does not anomalously mis-classify events by the 
sign of the charge. Performing the analysis without any 
selection on track charge in the FD produces consistent 
results. A scan by eye of events in the ND showed the 
track reconstruction efficiency to be well modeled. 

Differences in the underlying inelasticity distributions 
for and events lead to differences in the muon 
and hadron energy distributions for the two samples. 
Studies to validate the reconstruction of muon tracks 
and hadronic showers included the tightening of recon- 
struction quality criteria, investigation of the compara- 
tive performance of various detector regions, and the use 
of an alternative hadronic energy estimator. These stud- 
ies show the detectors to be well modeled, and that any 
mismodeling in reconstruction and selection efficiencies is 
accounted for in the systematic uncertainty. In conclu- 
sion, no evidence is found for any systematic error that 
could cause a significant difference between the measured 
i'^ and Vp oscillation parameters. 

In summary, MINOS has used a beam optimized for 
the production of 77^ to make the first direct observa- 
tion of disappearance and, assuming the disappear- 
ance is caused by oscillation, has made the most precise 
measurement of the corresponding antineutrino mass- 
squared difference to date. From fits to these data. 



MINOS has measured the oscillation parameters to be 
|Ato2| [3.36l^;4j^(stat) ± 0.06(syst)] x IQ-^eV^ and 
sm'^{20) ^ 0.86lE;;{^(stat) ±0.01(syst). The MINOS Vf, 
and measurements are consistent at the 2.0% confi- 
dence level, assuming identical underlying oscillation pa- 
rameters. Additional data are currently being taken with 
the F^-enhanced NuMI beam, with the aim of doubling 
the statistics presented in this Letter. 
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